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Abstract 

We investigate a particular type of curvaton mechanism, under which inflation can oc- 
cur at Hubble scale of order 1 TeV. The curvaton is a pseudo Nambu-Goldstone boson, 
whose order parameter increases after a phase transition during inflation, triggered by 
the gradual decrease of the Hubble scale. The mechanism is studied in the context of 
modular inflation, where the inflaton is a string axion. We show that the mechanism 
is successful for natural values of the model parameters, provided the phase transition 
occurs much earlier than the time when the cosmological scales exit the horizon. Also, 
it turns our that the radial mode for our curvaton must be a flaton field. 



Inflation is the only compelling theory to date for the solution of the horizon and flatness 
problems of the big bang cosmology as well as for explaining structure formation in the Universe. 
Recent precise observations have confirmed the basic predictions of the inflationary paradigm 
by ascertaining the spatial flatness of the Universe and the approximate scale invariance of the 
density perturbations, which give rise to the anisotropy of the Cosmic Microwave Background 
Radiation (CMBR) and seed structure formation. These exciting developments have rendered the 
inflationary paradigm a necessary extension to the hot big bang standard cosmology. 

In the light of precision data, inflation model-building can be upgraded beyond the simple 
single-fleld stage of its early beginnings. Indeed, more complex and realistic models of inflation, 
with tighter connections to the theory, less fine tunning and enhanced predictability and falsifi- 
ability are now possible to construct, making use of the rich content of particle physics. A first 
such example is the well-known hybrid inflation model [1], which couples the inflaton field to the 
Higgs field of a Grand Unified Theory (GUT) in order to obtain without tunning the desired false 
vacuum energy scale [2]. In hybrid infiation the inflationary period is terminated through the 
dynamics of this other fleld. 

In an analogous manner, one can attribute the generation of density perturbations during 
inflation to a field other than the infiaton [3]. This so-called curvaton field allows inflation to 
take place at a much lower energy scale than the typically required GUT-scale [4] and, in general, 
may relax a number of constraints regarding inflation model- building [5]. Low-scale inflation can 
revamp a number of inflation models that are well motivated on particle physics grounds [4]. It is 
important to stress here that the curvaton is not an ad hoc additional degree of freedom introduced 
"by hand" , but it may be a realistic field, already present in simple extensions of the standard 
model. Indeed, many such examples exist in the literature [6, 7]. 

However, even when a curvaton field is considered, there exists a lower bound for the infla- 
tionary scale, which, for generic curvaton models, can be quite tight [8]. This lower bound can 
be substantially relaxed for certain types of curvaton models [9], which enables inflation to be 
directly connected to realistic, beyond the standard model physics. 

In this letter I present a curvaton model which allows inflation at a Hubble scale as low as 
1 TeV. The curvaton field is a pseudo Nambu-Goldstone boson (PNGB), whose order parameter 
is substantially increased after the cosmological scales exit the horizon during inflation. As shown 
in [9], the result of this increase is to amplify the curvaton's perturbations. This enables even 
low-scale infiation to generate density perturbations of the observed amplitude. In the curvaton 
model presented, the increase of the PNGB order parameter follows a phase transition during 
infiation, which releases the radial mode from the top of the potential hill. 

The use of a PNGB curvaton is highly motivated because such a curvaton can be naturally 
light during inflation, since its mass is protected by the global U(l) symmetry [7]. This dispenses 
with the danger imposed by supcrgravity corrections, which typically lift the flatness of the scalar 
potential [10]. We investigate the performance of the curvaton model in the context of modular 
inflation, which corresponds to Hubble scale of order 1 TeV. Modular inflation is a well motivated 
model, which uses a string axion as the inflaton [11]. 

Let us begin by presenting the amplification mechanism for the curvature perturbations. We 
discuss here the case of an PNGB curvaton, whose order parameter has a different (larger) expec- 
tation value in the vacuum than during inflation and, in particular, when the cosmological scales 
exit the horizon. Thus, the potential for the curvaton field a is 

V{a) = {vm^f[l-cos{a/v)] V {\a\ < v) ^ ^mla^ , (1) 

where v = v{t) is the order parameter determined by the expectation value of the radial field \4>\ 
and rha = rha-iv) is the mass of the curvaton at a given moment. In the true vacuum we have 
V = vo and rha = rria with ?;o being the vacuum expectation value (VEV) of the radial field and 
nia being the mass of the curvaton in the vacuum. 

Let us demonstrate that the curvaton perturbations can be amplified by the non-trivial evolu- 
tion of the radial field. We begin by using the fact that [3] : 

C ^ ^decCa , (2) 
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where ( ~ = 2 x 10~^ is the curvature perturbation of the Universe, fidcc < 1 is the density 
fraction of the curvaton density over the density of the Universe at the time of the decay of the 
curvaton: 

Pa 



^dec 

p 



<1 (3) 

dec 

and C<7 is the curvature perturbation of the curvaton field a, which is given by 



Sa 

Ca 

a 



(4) 

dec 

where 'osc' denotes the time when the curvaton oscillations begin. Note that, non-gaussianity 
constraints from the observations from the WMAP satellite [12] restrict the range of il^ec 
follows: 

10-2 < Odec < 1 • (5) 

In this paper we consider the inflationary Hubble scale to be comparable to the tachyonic mass 
of the radial field, which determines the value of the order parameter of our PNGB curvaton. This 
means that the evolution of the radial field ceases at (or soon after) the end of inflation. Therefore, 
at the end of inflation, v ^ vq and the mass of the curvaton assumes its vacuum value mo-. Hence, 
in the following we assume that the curvaton mass has already assumed its vacuum value before 
the onset of the curvaton oscillations. Consequently, the curvaton oscillations begin when 

^^osc ~ ■ (6) 

Before the oscillations begin the curvaton is overdamped and remains frozen. This means that 
6*080 — where the denotes the values of quantities at the time when the cosmological scales 
exit the horizon during inflation and 

e = a/v, (7) 
with 9 e (— TT, tt]. Hence, for the curvaton fractional perturbation we find 

- = - :^ - . (8) 

<^ * ^ * osc 

Now, for the perturbation of the curvaton we have 

Sa. = ^, (9) 

We assume that the expectation value of the radial field during infiation is smaller compared to 
its VEV by a factor 

£ = - « 1 . (10) 

Combining Eqs. (8), (9) and (10), in view also of Eq. (7), we find 

5aosc ^ ^ , (11) 

which means that after the end of inflation, when the radial field assumes its VEV, the curvaton 
perturbation is amplified by a factor (see Figure 1). Prom Eqs. (2) and (4) we have 

CTosc ^ (^dec/O'^O-Qsc • (12) 

Using Eqs. (9) and (11), we can recast the above as 
We may obtain a lower bound on e as follows: 



(13) 



< 1 e> ffmin = 7, ' (14) 



where we have used Eqs. (7), (9) and (10) and that CTosc ^ vo- 
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Figure 1: Schematic representation of the amplification of the PNGB curvaton perturbation, when 
the order parameter v increases from the vahie it has when the cosmological scales exit the horizon 
= evo to its vacuum value vq- The perturbation at horizon crossing has amplitude 5(T* ~ i?*, 
which corresponds to a phase perturbation for the radial field |(/)| of magnitude 69 = Sa^/v^. As 
the order parameter grows 69 remains constant (the phase perturbation is frozen on superhorizon 
scales) but the amplitude of the curvaton perturbation is increased up to Sa ~ e^^i/*. 



Now, as is shown in [9], in the case when the curvaton oscillations begin after the radial field 
has attained its VEV, we have"'^ 

- f^Zef ( — rr^V^V' ( "'"^!;""''^-^ V''(vreC)^/^(r|3.r.3.)V5^ (15) 



min{TO<^,ri„f}/ V Hj 

where Finf and F^ are the decay rates of the inflaton and the curvaton fields respectively, -ffdom 
is the Hubble parameter at the time when the curvaton density dominates the Universe (if the 
curvaton does not decay earlier) and i?BBN ^ T'^B^I^np is the Hubble parameter at the time of Big 
Bang Nucleosynthesis (BBN), with Tbbn ^ 1 MeV. 

Now, we require that the curvaton decays before BBN, i.e. F^. > H^^^. We also have Finf < H^. 
Hence, Eq. (15) provides the following bound 

H* > ^ioJ'i^^O^^HTl,m%y^' ^ (eV^^dcc)'/' x lO^GeV . (16) 
Furthermore, we also note that 

r.>^, (17) 



2 

mi 



where the equality corresponds to gravitational decay. The above can be shown [9] to imply that 

> ^-.LineCfmp .^axjl,^!'^', (18) 



^We use natural units, where c = h = 1 and Newton's gravitational constant is G = dnrmp^ with 
mp = 2.44 X lO^^GeV being the reduced Planck mass. 
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which results in the bound 



> n^,\{neCfrnp (^^^ ~ (eV^^dec) x lO^OCeV . (19) 

This bound may be relaxed if e is small enough. In particular, for a PNGB curvaton we may have 
< if*. Comparing the bound in Eq. (16) with the one in Eq. (19) we find that the former 
bound is more stringent if 



^ < (—] (—) ~ 10-^ n-ll\H./m.f/' (20) 

Thus, for e <C 1, the second bound is typically less stringent than the first one. 

During inflation, the evolution of the order parameter of the PNGB curvaton, is subject to an 
important constraint, which has to do with preserving the scale invariance of the spectrum of the 
curvature perturbations. 

The amplitude of the density perturbations is determined by the magnitude of the perturba- 
tions of the curvaton field, which, in this scenario, apart from the scale of if* is also determined by 
the amplification factor e~^. The latter is determined by the value of the order parameter when 
the curvaton quantum fluctuations exit the horizon during inflation. A strong variation of v{t) at 
that time results in a strong dependence of e(fc) on the comoving momentum scale k, which would 
reflect itself on the perturbation spectrum threatening significant departure from scale invariance. 

In Rcf. [9] is was shown that, in order for this to be avoided, the rate of change of the radial 
field must be constrained as 

\v/v\^ = l^/^l* < , (21) 

where |0| is the radial field, which determines the value of the order parameter. In fact, the 

contribution to the spectral index due to the evolution of v is Sus = —2H~^{v /v)^,. From the 
above it is evident that, in order not to violate the observational constraints regarding the scale 
invariance of the density perturbation spectrum, the roll of the radial field has to he at most very 
slow when the cosmological scales exit the horizon. However, this cannot remain so indefinitely 
because we need vq ^ v^, to have substantial amplification of the perturbations (i.e. e <C 1). 
Consequently, v has to increase dramatically at some point after the exit of the cosmological 
scales from the horizon. This requirement is crucial for model-building.^ 

In our model, we will show that the evolution of v begins at a phase transition during inflation. 
Initially, the growth of v is very slow, but later, near the end of inflation, v grows substantially 
until it reaches its vacuum value vq. 

Let us now briefly describe the model of inflation. We are going to consider a PNGB curvaton 
a whose radial field |0| is of bare mass similar to the Hubble parameter during inflation, that is 

~ if* . (22) 

This has the advantage that the radial field rolls substantially by the end of inflation so that s can 
be very small. In particular, we will assume that the tachyonic mass of the radial field is a soft 
mass generated by supersymmetry breaking and it is, therefore, roughly of the electroweak scale 
7713/2- Hence, we consider infiation at the intermediate scale 

y//* ~ V^3/2mp ~ IQio-^GeV if* ~ rns/s . (23) 

A particular example of such an inflation model (but, by all means, not the only one) is modular 
inflation [11], where the inflaton field s is a string axion, whose flatness is lifted by gravity mediated 
supersymmetry breaking. In this model the inflationary potential is of the form: 

V{s) = Vinf - + ■■■ , (24) 



^The requirement in Eq. (21) may be even more fundamental in origin. Indeed, a PNGB with rapidly varying 
order parameter cannot be treated as an effectively free field. I would like to thank D.H. Lyth for pointing this out. 
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where the elhpsis denotes terms, which are expected to stabiUse the potential at Svev ^ Trip. 
Therefore, in the above we have 

Vinf {m3/2mpf and Hi^f 7713/2 , (25) 

where Hmf — VMnf /3mp. 

This inflation model results in fast roll inflation [13], where 

s = Si„exp(^^,AiV) and F^-^ ( ^1 + 4c/9 - l) with c= (26) 

where AiV is the number of the elapsed c-foldings. From the above one can easily obtain the 
inflation scale e-foldings before the end of inflation, as 

V{N) ~ Vinf (1 - e-^^^^) . (27) 

Even though fast-roll, modular inflation keeps the Hubble parameter H rather rigid. Indeed, 

it can be easily shown that 

because ~ 1 and s <^ nip during inflation, with e <C 1 being one of the so-called slow roll 
parameters defined as 

e^--J. (29) 

For modular inflation the initial conditions for the inflaton fleld are determined by the quantum 
fluctuations, which send the field off the top of the potential hill. (The modulus can be considered 
to be originally placed at the local maximum because the latter can be thought of as a fixed point 
of the symmetries.) Hence, we expect that the initial value for the inflaton is 

Si„~i?inf/27r. (30) 

Using the above and considering that the flnal value is Svev ~ ^np, we can estimate, through the 
use of Eq. (26), the total number of e-foldings as 

A^tot-^mf^V (31) 



"T-3/2 



where we took into account Eq. (25). 

Let us turn our attention to the curvaton model. Consider the superpotential^ 

W=^^^, (32) 
n + 3 ^ ^ 

where n >0 and the complex field (j) can be thought to contain the curvaton phase field a and 
one radial field as follows: 

(j) = \(j)\e^^ = |<^| exp{a/V2v) . (33) 
Then the scalar potential can be written as 



A 



{CaH + A) " „ +/I.C. 



2n+4 



{C,H- - m^)|0p + A2^L_ + ^CaH + ^)^M^ cos[(n + 3)0] (34) 



^Such type of superpotential is reminiscent of supersymmetric realisations of the Peccei-Quinn symmetry in 
which the Peccei-Quinn scale is generated dynamically [14]. 
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where and A are soft supersymnictry breaking mass-scales at zero temperature, both given 
by the electroweak scale m3/2- Note that we have put negative mass-squared for the field 
at zero temperature to break the U(l) symmetry. We also considered corrections coming from 
supergravity, which provide effective mass terms of order H [10] (for their effect on curvaton 
physics see Ref. [15]). Absorbing the (n + 3) factor into 9 (and shifting the latter by tt) ^ we can 
write the curvaton potential as: 



V{a)K\{CAH + A)v^ (^-^^ [l_cos(0] . (35) 



We are going to assume that the U(f ) symmetry is broken at some moment during inflation 
with ~ TO0 ^ ^3/2- Hence we take ^ +1. After this moment the radial field begins to 
grow, which can result in £ <C 1. In time, after the symmetry breaking, the tachyonic effective 
mass of the radial field approaches its vacuum value as the supergravity correction diminishes 
due to the gradual decrease of the Hubble parameter. 

After the phase transition, the time-dependent minimum of the potential of the radial field is 
given by 



|.^|min = (A-im^y'ml - C4,m) , (36) 

which gradually grows. Soon l^jmin assumes its vacuum value: 

t-o- (A-V^m^)^, (37) 
From the above and also in view of Eqs. (1) and (35) we find 

ml^X{CAH + A)v(—\ . (38) 



Evaluating the above after the order parameter assumes its vacuum value v vq we obtain 

ml « {CaH + A)m^ , (39) 
where we used Eq. (37). Since ^ 7713/2, Ca ^ 1 and H < H^, we find that 

rria ^7713/2^^/"*. (40) 

However, during inflation the effective mass of the curvaton is much smaller. Indeed, in view 
of Eq. (38), we get 

— I rhcriv*) e~ rria , (41) 

voj 

where we used Eq. (10). Therefore, since £ <^ 1 and ^ 7713/2 ^ H^, we see that, during inflation 
7710- «C -ff*, i.e. the PNGB is appropriately light and can act as a curvaton field. 

Let us now calculate the value of e required so that the scenario works. Firstly, we note that, in 
our case, the curvaton assumes a random value at the phase transition, which typically is a <^ v. 
After the end of inflation and before the onset of the oscillations the field is overdamped and 
remains frozen. Hence, we expect that at the onset of the oscillations we have: 

CTosc Ovq , (42) 

where, typically, 9 ^ 0(1) and we took into accoTint that the radial field assumes its VEV very 
soon after the end of inflation. Combining Eqs. (13) and (42), we find 

-%(^)*, (43) 
Tr(9 \ nip J 



''This, in effect, means considering the range: < — n < 
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where we also used Eq. (37) taking m3/2 and A 1. The above is always larger than £ 

where 

' m3/2 \ ^ 
mp J 

where wc also used Eq. (14) with ^ rn^^2- 

Let us now enforce the constraint in Eq. (16), which, for if* ~ reads 



(44) 



e<^(^) (^] ~ lO-v/lW. (45) 



ttC VTbbn/ V mp 

Prom Eqs. (43) and (45) it is easy to find that the above bound can be satisfied only if n is large 
enough: 

^ 8 + iog(v/rw/^) 

''>7-iog(yrw/0)- ^'"^ 

According to Eq. (5), we see that, at the best of cases, (when Ojec ^ 10~^ and ~ 1) we have 
n > 1. Hence, we see that the radial field must correspond to a flaton field, stabilised by non- 
renormalisable terms. An upper bound on n can be obtained by requiring that the curvaton 
decays before BBN. 

The decay of the curvaton depends on its coupling to other particles. The lowest decay rate 
corresponds to gravitational decay with ^ m^/m^p. However, if 4> is part of a supersymmet- 
ric theory we may expect a much larger value for An interesting possibility is realised by 
introducing the following coupling between (j) and the Higgses: 

In this case, as is evident from Eq. (37), our curvaton model also solves the /U-problem for A/j/A ~ 1. 

Now, the interaction of a with ordinary particles is governed by the effective /i-term in Eq. (47), 
which results into the following decay rate of a into two Higgs particles: 

(«^^ 

Demanding that > H^bk results in the bound 

~ 10-^ TeV > i^BBN ~ lO-^^TeV => m„> lO^TeV , (49) 

where we used Eq. (37). Since we consider m^ ~ mz/2 < TeV we see that there is a mild upper 

bound on n which, roughly, demands n < 9. 

To proceed further, we have to consider separately the cases when the curvaton decays before 
or after it dominates the Universe. Suppose, at first, that the curvaton decays before domination 
(fldec ^1)- In this case, during the radiation epoch and after the onset of the oscillations, for the 
curvaton density fraction we have Pa/ P oc a{t) cx iJ^^/^. Hence, we find 

• r 1 X 1/2 / 

mm{TOCT,l infj \ <J. 



^^dec- """^--'-""^ — , (50) 



where we have used Eq. (13) and 



2 

Pa_ 
P 



mp 



(51) 



with {pa)osc — ^"^CT^'osc ^^'^ Pose ^ m^m^,. Using Eq. (48) into Eq. (50) and also Eqs. (37) and 



2 

(13) we obtain 



(52) 

ttC V mp J 
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where we have also used that Finf < 1^3/2 ^ n^a- and 



Tinf ~ g 1713/2 , (53) 

with the mass of the inflaton field s taken to be < ff* ~ W3/2 and g being the couphng of 
the inflaton to its decay products. In principle, g can be as low as nis/mp if the inflaton decays 
gravitationally. However, since reheating has to occur before BBN, g has to lie in the range: 

10-14^ 10 !^^<5<1. (54) 
Combining Eqs. (43) and (52) we find the relation 

' ("^"^ ^ (55) 



Let us now consider the case when the curvaton decays after domination (Odoc ~ 1). In this 
case, the curvaton dominates the energy density of the Universe when H = .ffdom, where .ffdom is 
given by 

^^dom ^ ( ) min{m„,ri„f} . (56) 
\mp J 

Now, using Eqs. (42), (48) and (53) it can be shown that the requirement r,^ < .ffdom results in 
the bound 

The fact that the case of curvaton domination requires a larger value of g [compare the above with 
Eq. (55)] is to be expected because, this means that the inflaton decays earlier and, therefore, the 
density fraction p„/ p grows substantially, allowing the latter to dominate the Universe before its 
decay. The higher g is the more dominant the curvaton will be. 

Note also, that, when the curvaton decays after it dominates the Universe, the hot big bang 
begins after curvaton decay, which suggests the reheating temperature 

T'reh ~ V^amp 1713/2 I 1 . (58) 

It can be easily checked that the above is higher that Tqbn when n < 9, in agreement with the 
bound from Eq. (49). 

From Eqs. (55) and (57) we see that, in general. 

For ^ 1 and in view of Eqs. (5) and (54) the above bound suggests 

n > 2 , (60) 

which is tighter than the bound in Eq. (46). 

We now concentrate of the evolution of the radial field \(f)\, which has to be such, as to achieve 

the required value for e. Let us assume, at first, that the radial field follows the growth of the 
temporal minimum given by Eq. (36). In this case we can calculate the amplification factor as 



1)* 



(61) 
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Using Eqs. (38) and (36) one finds that, in this case, the curvaton's mass is given by 

ml « {CaH + A) - C^H^ (62) 

The above and (61) suggest that 

{^^\[A + CAH.)^e-+\ (63) 

which agrees with Eq. (41), given that Ca ~ 1 and i/* ~ A ~ ~ ma ~ m^/2- From Eq. (36) 
it is easy to show that the rate of growth of the order parameter is 

V |</)U„ n+l[c^H^ ') ^^^^ 

where e is defined in Eq. (28). From Eqs. (61) and (64) we obtain 

(u/t;), (65) 

Comparing this with Eq. (21), we find that, for the scale invariance of the spectrum to be preserved, 
we require 

<£2(n+i)^ (66) 

where = e(s*). 

Now, if the growth of |0|min is so rapid that the radial field cannot follow it, then we expect 
to roll, instead, down the potential hill. In this case the order parameter is determined by the 
rolling 

When the cosmological scales exit the horizon the radial field has to be slowly rolling because 
we need the order parameter to vary slowly enough, not to destabilise the approximate scale 
invariance of the perturbation spectrum [cf. Eq. (21)]. Therefore, the Klein-Gordon equation for 
m is: 

3H,\<P\-fnl\4>\c^0, (67) 

where 

ml = ml-C4,H\ (68) 
Using the above, the rate of growth of the order parameter, in this case, can be easily found to be 

The variation of the order parameter is expected to follow the less rapidly changing rate of 
growth. Hence, by comparing the two rates in Eqs. (64) and (69), we see that the order parameter 
follows the variation of |(?i|min only if 

> e* , (70) 

where we used again Eq. (61) and also ~ 1. It is evident that, if the above constraint is satisfied 
then so is the requirement in Eq. (66). Note, however, that if the above constraint is violated then 

the order parameter v is determined by the rolling |(/)| and not by the varying |(/>|min, in which case 
the requirement in Eq. (66) is not valid, while also the amplification factor is not the one shown 
in Eq. (61).^ 

In this latter case, we find the amplification factor as follows. Using Eq. (27) we can write \cf)\ 
as a function of the number A'^ of the remaining e-foldings of inflation. Starting from Eq. (67) and 
after a little algebra we obtain 

3 rflnl^l _ e-^''^^^ - c-'^'-^ 



Crf, dN 1 



^If V follows the growth of \(j>\ instead of |</i|niin then e is expected to be smaller that the one in Eq. (61) because, 
at any given time, v{t) < |</>|min(t)- 



9 



where ./Vx corresponds to the phase transition which changes the sign of m|. By definition 

ml ^ C^Hl C^HUI - e-2^»^>-). (72) 
where = H{N^). Integrating Eq. (71) we get 



- 2{N^ - N,), (73) 



where |0|x = MN^)- 

The displacement of the field from the origin at the phase transition is determined by its 
quantum fluctuations. This means that 

\ct>\^:^H^/27T, (74) 

We also have 

s=\^ = \^^^\^U^^^, (75) 

Vo -H* Vo £rnin StT 

where we used Eq. (14). 

In view of Eq. (69), the requirement in Eq. (21) becomes 

where wc took into account Eq. (72). 

Finally, another issue to be addressed concerns the requirement that the radial field does slow 
roll at the time when the cosmological scales exit the horizon. In order for this to occur, its 
quantum fluctuations should not dominate its motion, i.e. |(/)| has to be outside the quantum 
diffusion zone. The condition for this to occur is ff*/27r < {(p/H)^, or equivalently 



dV 



2fhl{H.)\<PU > (77) 



Using Eqs. (68) and (72) and working as before, the above constraint is recast as 

(m) ~ '"^^'/"^ ^ (t^^^SSvo ) + 1 {t^^S:) ■ 

To illustrate the above we present an example, taking 

n = 2 and A, 6* ~ 1 . (79) 

The bound in Eq. (49) suggests that this is acceptable provided rua- > 5 GeV. Using Eq. (43) we 
obtain the value of the amplification factor, necessary for the model to work: 

s ~ 10-^Odec . (80) 

If the curvaton decays after domination then Eq. (57) demands .g > 1, which is not compatible 
with the range in Eq. (54). Therefore, we have to assume that the curvaton decays before domina- 
tion, in which case fidec ^ 1, with the bound saturated when the curvaton decays approximately 
when it is about to dominate the Universe. In this case, Eq. (55) suggests 

g ~ fidec < 1 . (81) 

Such a large coupling can be understood only if the VEV of the inflaton modulus is an enhanced 
symmetry point. As a result of the above, the reheating temperature after the end of inflation is 
found to be 

Treh ~ 5(V™3/2mp ~ fi( X 10^°'^ GcV. (82) 
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From the above we see that, in order not to challenge the gravitino constraint, we have to choose 
the lowest possible value of g, which, according to Eqs. (5) and (81) corresponds to 

fidec ~ 10-^. (83) 

Hence, from Eqs. (80) and (81) we obtain the values 

e ~ 10"^ and g ~ 10"^. (84) 

From Eq. (30) and (25) and also, using Eq. (28), it is easy to see that 

e* > 10-30, (85) 

where wc considered that > Sjn. Hence, from Eqs. (84) and (85) it is straightforward to sec that 
the constraint in Eq. (70) is badly violated, which means that the order parameter v follows the 
slow roll of the \(f)\ field and not the variation of the minimum of the potential |(/)|min- Consequently, 
the amplification factor e is not given by the; expression in Eq. (61) in this case. Instead, we can 
estimate the amplification factor with the use of Eq. (75). 
Using Eq. (44) with n = 2 we find 

- 10-1°. (86) 

Therefore, Eqs. (44), (75) and (84) suggest 

l</'U~10'§. (87) 

The above can, in principle, be used in Eqs. (73) and (78) to constrain the parameters of the 
underlying model. 

A useful quantity to calculate in order to evaluate Eqs. (73) and (78) is the number of e-foldings, 
which corresponds to the cosmological scales iV* . The cosmological scales range from a few times 
the size of the horizon today ~ H^^ down to scales ~ 1Q~^Hq^ corresponding to masses of order 
IO^Mq [16]. Typically this spans about 13 e-foldings of inflation. For the estimate of iV* we will 
chose a scale roughly in the middle of this range; the scale that re-enters the horizon at the time 
when structure formation begins, i.e. at the time teq of matter-radiation equality. Then, in the 
case when the curvaton decays before domination it is straightforward to obtain 

,1/3^1/6 



exp(7V0 ~ Ht''T{i°^^^ fl'^'V^^^^M^ (88) 

where we have used Eq. (53) and that ff* ~ mo- ~ ^3/2- Using Eq. (84), we obtain 

iV* ~ 43 . (89) 

The number of c-folds that corresponds to decoupling (when the CMBR is emitted) is roughly 
+ 1.5, while the one which corresponds to the present horizon is ~ A^* -|- 9. 
In the attempt to obtain the allowed parameter space for our model it soon becomes clear that, 
while the requirement in Eq. (78) is relatively easy to satisfy, the major difficulty is reconciling 
Eq. (73) with the bound in Eq. (76) coming from the spectral index requirements. This is especially 
true in view of the recent WMAP results [17], which correspond to spectral index = 0.96 ± 0.02, 
i.e. Us > 0.92 at 95% c.l. This means that the left-hand-side of Eq. (76) should not exceed 0.04.^ 
By careful investigation of Eqs. (73) and (76) it is found that the above difficulty is more alleviated 
the larger the value of A'x is, i.e. the earlier the phase transition occurs. In fact, a solution is only 
possible if 

2F,N^ » 1 . (90) 

®Note that, in our model, all other contributions to the deviation of the spectral index from unity [18] are 
negligible. 
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In view of the above Eqs. (73) and (76) can be respectively approximated as 



< 0.12 (e^^^^* -1). (92) 

Now, using Eq. (74) we can write: 

where we have considered also Eq. (27), using that H^{N) ~ V{N)/3mp. In view of the above 
and according to the approximation in Eq. (90) we can recast Eq. (91) as 

mm^-rU^^Vnd-e— ). (94) 



(95) 



Under the same approximation Eq. (78) becomes 



In (^^IlMl^ > 2F,N. - HC^M + ln(l - e-'^'^'] 



where we have also used Eq. (93). 

Solving Eq. (94) in terms of and using Eq. (92) we obtain 



2£yv« 

86 



_^ 4 In 10 



ln(l - 



2 F.N,-' 



0.04(e2^»^* - 1) < (96) 



where we have also employed Eqs. (87) and (89). Solving numerically we obtain the bound 

^'a < ^ - 1-8 X 10-3 (97) 
This bound, in view of Eq. (26) results in 

< 0.073 i/inf . (98) 

which is somewhat tight and implies that inflation is not really of the fast-roll typo, but the inflaton 
is light enough to roll slowly down its potential hill. From Eqs. (31) and (97) one obtains 

iVtot > 1-9 X 10^ (99) 

Thus, if the phase transition, which releases |(/)| from the origin, occurs not much later than 
the onset of inflation, then the approximation in Eq. (90) can be well justified. Similarly, using 
Eqs. (89) and (97), Eq. (92) gives the bound 

< 0.020 . (100) 

In view of Eqs. (72) and (90) the above bound suggests 

m^<0.14iri„f. (101) 

The values of rus and can approach iJinf if one decreases but then the constraint in Eq. (76) 
becomes seriously challenged. It can be easily checked that, with the above values the requirement 
in Eq. (95) is satisfied as well. 

The above results suggest that, for the n = 2 case and when A, ^ 1, the model can work for 
masses of the order 

m^<0.1m3/2 and mg < 0.01 m3/2 , (102) 
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where 77x3/2 ~ 1 TeV. Such values imply only a mild tuning on the masses: predominantly on the 
mass of the inflaton modulus. This is necessary because the variation of H should be kept small, 
since only then can the tachyonic effective mass of the radial field remain small enough for \(f)\ 
to be slow-rolling and the constraint in Eq. (21) to be satisfied. Note that, a tuning of the inflaton 
mass is quite plausible, since the latter is a string axion. 

One may wonder why, since both the inflaton field s and the radial field turn-out to be light 
when the cosmological scales exit the horizon during inflation, we cannot use those fields to generate 
the observed curvature perturbations. The reason is that, in contrast to the PNGB curvaton, the 
perturbations of those fields are not amplified. Hence their contribution to the overall curvature 
perturbation is insignificant. Indeed, for the inflaton we have ^ (nis/s^,) ^ 10^"'^^, which is 
much smaller than the observed value C — 2 x 10~^. Similarly, for \<p\ it is easy to show that 
C<p ^ sCcr ^ 10~^^, where we used that w (. 

In conclusion we have seen that our mechanism can work with natural values of the parameters 
with only a mild tuning on the inflaton mass. Another important requirement is that the phase 
transition, which releases the radial field from the origin, occurs much earlier than the time when 
the cosmological scales exit the horizon, in order not to destabilise the flatness of the curvature 
perturbation spectrum. 

Our PNGB curvaton is such that can be easily accommodated in simple extensions of the 
standard model. Indeed, in Ref. [19] we present in detail such a realisation, using as curvaton an 
angular degree of freedom orthogonal to the QCD axion in a class of supersymmetric constructions 
of the Peccei-Quinn symmetry. Presumably, other PNGB curvatons, such as the ones in Ref. [7], 
can also be utilised. 

We should note here that, although the modular inflation model, which we considered, is highly 
motivated, it is by no means the only possibility. Other inflationary models with Hubble-scale of 
order 1 TeV may also be applied [20] . Needless to say that designing inflationary models at such 
energy scale can allow direct contact with particle physics. 

I am grateful to G. Lazarides for stimulating discussions. 
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